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Abstract. 



The suppression of hadron p T spectra in high energy central heavy-ion collisions com- 
pared to proton-proton collisions, referred to as 'jet-quenching', is currently attributed to 
partonic energy loss in the hot medium created in the collision. The RHIC experiments show 
that at large enough p T , hadron quenching is strong, and of similar magnitude for light and 
heavy flavours. This point is difficult to understand in a parton energy loss scenario, where 
the energy loss is believed to be dominantly radiative, and quantitatively different for light 
partons and heavy quarks: gluon radiation off a heavy quark is suppressed at small angles - 
within the so-called dead cone - and a heavy quark suffers less radiative energy loss than a 
light parton. In this context it is important to reconsider the collisional contribution to par- 
tonic energy loss. Although it seems difficult to see how collisional losses could substantially 
increase heavy flavour quenching without simultaneously increasing light hadron quench- 
ing, it is theoretically important to establish correct results for the heavy quark collisional 
energy loss. 
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At the Relativistic Heavy Ion Collider (RHIC), the p T spectrum of light hadrons 
has been observed [1, 2] to be strongly suppressed in central ultrarelativistic heavy- 
ion collisions, as compared to proton-proton collisions. This is illustrated by the 
factor R AA (p T ), 



appropriately normalized by the number iV co ii of binary nucleon-nucleon collisions 
in a nucleus-nucleus collision, to be much smaller than unity, as can be seen on 
Fig. 1. For pions R\ A ~ 0.2 for p T > 4GeV. This nuclear attenuation of hadron 
production, R h AA < 1, has been quoted as jet- quenching, and is one of the most 
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FIGURE 1. Light hadron quenching. Direct photons are not quenched {FC AA ~ 1), whereas 
light hadrons are (i?^ - 0.2). Figure taken from Ref. [3]. 



important signals for the new state of matter - referred to below as a quark-gluon 
plasma or QGP - produced in a (central) heavy-ion collision. 

The effect of jet-quenching was anticipated by Bjorken [4], as a consequence of 
medium-induced partonic energy loss in the QGP. Although Bjorken illustrated 
his idea by performing an estimate of the collisional loss of light partons, parton 
energy loss generally arises from both collisional and radiative processes. In fact, 
for a parton energy E ^> M,T, where M is the parton mass and T the plasma 
temperature, the radiative component is believed to be dominant. For a medium 
of large size we have [5] 

AE rad (L > L cr ) ~ a]^ET^L > AE cM ~ a 2 s T 2 L . (2) 
The critical length L cr reads [5] 

1 / E 1 

L cr -W->A -, (3) 

a s T V T a s T 

where A is the parton's mean free path between two elastic scatterings. When L is 
not large, L < L cr , the radiative loss is modified [6, 7], 

AE rad (L < L cr ) ~ a 3 s T 3 L 2 , (4) 

but still dominates over the collisional loss as long as L > X. Thus it seems 
legitimate to neglect collisional energy loss in phenomenological studies of light 
hadron quenching, and to attribute hadron quenching solely to the radiative parton 
energy loss 1 . 



For a light quark and for L = 5fm, a careful comparison of collisional and radiative contributions 
within the same consistent model [8] gives the estimate AE co u / AE ra d ~ 0.25 at E ~ 20 GeV. Using 
A ~ lfm, this is in agreement with our estimate AE co u / AE rac i ~ X/L when L < L cr . In Ref. [9] 
a surprisingly large value AE co u / AE rac i ~ 0.75 is used at E ~ 20 GeV, obviously changing the 
qualitative understanding of quenching. 
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FIGURE 2. The Raa factor for electrons from heavy D and B meson semi-leptonic decays. 
Figure taken from Ref. [12]. 



The data on heavy flavour quenching [10, 11], obtained indirectly by observing 
the electrons and positrons arising from D and B meson semi-leptonic decays, shows 
that the nuclear attenuation of heavy c and b quarks is as strong as that of light 
partons, Raa(Q) — RAA^.g) ~ 0.2 (see Fig. 2). This is difficult to explain when 
assuming that heavy quarks, similarly to light partons, lose energy dominantly 
through radiation. Indeed, for heavy quarks (M 3> Aq CD ), gluon radiation at an 
angle 9 rad < M/pt is suppressed [13]. This is the so-called 'dead cone effect', leading 
to a strong suppression (at not too large pr), of the average radiative loss of a heavy 
quark compared to a light parton, 

AE rad (b) < AE rad (c) < AE rad {q, g) . (5) 

If the electrons and positrons from heavy flavour decays observed in the experiment 
would arise dominantly from D mesons, then the dead cone effect would not be too 
drastic, due the moderate value of the charm quark mass. In this case explaining 
the heavy flavour quenching with purely radiative c quark energy loss might still 
be possible [8] . However, the relative contributions from D and B decays are likely 
to be of the same order, as in pp collisions [14]. The dead cone effect being stronger 
for b than for c quarks, this leads to a reduction of the effective radiative loss. 
As a consequence, with similar contributions from c and b quarks, the current 
phenomenological models typically underestimate heavy flavour quenching. 

Thus for heavy quarks, purely radiative energy loss seems insufficient to explain 
the observed attenuation. This has renewed the interest in the collisional part 
AE co u of the parton energy loss. Could the collisional loss be much larger for a 
heavy quark than for a light parton, which would increase heavy flavour quenching 
without spoiling the reasonably good description of light hadron quenching with 
purely radiative loss? 

We have reconsidered the average collisional loss of a heavy quark, and of a light 
parton. In the next section we summarize our calculations, which include the effect 



of the running coupling. We will see that the very definition of parton energy loss 
is different in the tagged heavy quark case and in the untagged light parton case. 
As a result, in the high energy limit E 3> M 2 /T, the heavy quark collisional loss 
contains an additional term - a collinear logarithm oc log(ET/M 2 ) arising from 
■u-channel exchange - which had been previously overlooked. This term makes the 
heavy quark collisional loss (25) slightly larger than the light quark loss (26), but 
is however small compared to the well-known soft logarithm oc \og{ET/m 2 D ) (mo 
being the Debye mass of the QCD plasma) arising from t-channel scattering. This 
is mainly due to a different associated color factor. We can thus approximate (see 
also (27)) 

E > M 2 /T AE coU (Q) ~ AE coll (q) . (6) 

We stress that the latter estimate holds in the ultra-high energy limit E ^> M 2 /T. 
However, taking E ~ pr — 20 GeV for the heavy quark energy, and T ~ 200 MeV, 
we see that the data represented in Fig. 2 correspond to the region E ^> M 2 /T 
for the charm quark, but not for the bottom quark for which E <C M 2 /T. Since 
AEcoii(Q) is an increasing function of E we must have 

M 2 /T » E » M 2 /T AE cM (b) < AE cM (c) ~ AE coll {q) . (7) 

Thus in the kinematical region of Fig. 2, we expect both collisional and radiative 
contributions to the b quark energy loss to be suppressed compared to the c or light 
quark case. In this respect the observed strong heavy flavour quenching, with an 
important contribution from b quarks, remains a puzzle. 

This conclusion should however be considered as preliminary. Indeed, an average 
loss smaller for b than for light quarks is a priori not incompatible with heavy and 
light flavour quenching being of similar magnitudes, due to the different energy 
loss probability distributions in both cases. 

2. PARTON COLLISIONAL ENERGY LOSS 

The collisional energy loss rate dE/dx of an incoming parton of 4- momentum 
P = (E,p) due to scattering off thermal particles of type % (i — q,g) is given by 

where we denote J k = J d 3 k/(27i) 3 . The expression (8) represents the thermal 
average - with appropriate Bose-Einstein or Fermi-Dirac thermal distributions - 
of the elastic scattering amplitude squared |.Mj| 2 , weighted by the energy loss in 
a single elastic scattering AE (to be defined precisely below). We will consider a 
fast parton of velocity v — > 1. The appropriate sum and average of \Aii\ 2 over spin 
and color is implicit in (8) and in the following. 
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FIGURE 3. Amplitudes for heavy quark elastic scattering in a QGP. A curly line denotes a 
gluon. The blob in (a) and (d) denotes the resummed hard thermal loop gluon propagator, which 
is necessary to screen the t-channcl contribution in the infrared. 



2.1. Tagged heavy quark 

In the case the incoming particle is a heavy quark of mass M ^> T, we can neglect 
the presence of quarks and antiquarks of the same flavour in the plasma, and the 
scattering amplitude M.i is given by the diagrams shown in Fig. 3. In the following 
we will denote P' = (E',p'), K = (k,k), K' = (k',k'), and use the Mandelstam 
invariants s = (P + K) 2 , u = (P-K') 2 , andt = Q 2 , where Q — K' - K — (u,q). We 
will focus on the limit E > M 2 /T > M, implying s = M 2 + 2P.K ~ O (ET) > M 2 . 

When the heavy quark is tagged in the final state, the energy loss AE appearing 
in (8) is simply the difference between the heavy quark initial and final energies, 
AE = E — E' = u. We stress that for a tagged particle, events corresponding to full 
stopping, i.e., to E' <C E or equivalently to maximal energy transfer u ~ cj max ~ E, 
contribute to dE/dx. 

When E ^> T the expression (8) can be simplified to [15] 

dE , f „, /" |t|max „ . do-i \t\ 

where di is the degeneracy factor of target particles of type i, and 

do~i 1 1 2 

~dt = 16tt(s — M 2 ) 2 d~d~i ^ ' (10) 

is the corresponding differential cross section. In Eq. (9) the bounds on t are set 
by kinematics, |t| max = (s — M 2 ) 2 / s ~ s, and |t| m i n = 0. 

The result (9) is exact in the limit E^T and can actually be simply understood. 
Noting that \t\ = 2K.K' = 2K.Q = 2(kcu - k.q), the factor AE to be inserted in (8) 
reads AE = uo = \t\/(2k) + k.q/ k. From the result (9), we see that effectively we can 
write uj = \t\/(2k), the term oc q.k vanishing after the angular integrals performed 
to obtain (9) starting from (8). Without the last factor uj = \t\/(2k) in the r.h.s. of 
(9) we recognize pa ~ 1/A, with p the spatial density of thermal scatterers, a the 
total elastic scattering cross section, and A the incoming particle's mean free path. 
Thus we can read the r.h.s. of (9) as ~ AE/X ~ dE/dx. 

In the following we will discuss, in a heuristic way, the leading logarithms 
contributing to dE/dx. We refer to Refs. [15, 16] for technical details. 



soft logarithm from t-channel exchange 



When \t\ <C s the t-channel diagrams of Figs. 3a and 3d contribute to the squared 
amplitude as \Mi\ 2 oc a 2 s s 2 /t 2 . Using (10), we obtain from (9) the contribution 
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We used \t\ mgx ~ s ~ ET, which is justified to logarithmic accuracy, and |t| m i n ~ m o 
as an effective infrared cut-off, where m 2 D = 47ra s T 2 (l +n//6) is the Debye screening 
mass in the QGP [17]. The logarithm in (11) thus arises from the broad logarithmic 
domain m 2 D <C |t| <C s. 

The inequality \t \ <C s implies that the energy transfer in a single collision satisfies 
ou <^ E. Thus the logarithm in (11) arises from kinematical configurations where 
the heavy quark is the leading parton in the 2-parton final state of the elastic 
scattering represented in Fig. 3. 



collinear logarithm from w-channel exchange 



The contribution from the square of the -u-channel amplitude shown in Fig. 3c 
brings another logarithm in the heavy quark energy loss, arising from the domain 
M 2 -C \u\ <C s ~ ET. Indeed, in this region we can approximate 
and \t\ ~ s (recall that s + u + t ■ 
thus obtain the contribution 
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2M ). Changing variables from t to u in (9), we 
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(12) 



We emphasize that the logarithmic contribution (12) arises from \t\ ~ s, or 
iv ~ c^max — E — M ~ E. Hence, in contrast to the soft logarithm of (11), the collinear 
logarithm of (12) arises from full stopping of the incoming heavy quark, where the 
leading particle in the final state of Fig. 3 is the struck thermal parton, to which the 
incoming heavy quark has tranferred all its energy. This effect is the QCD analogue 
of QED Compton backscattering of laser beams, used to produce energetic photons 
from energetic electrons. Full stopping must be taken into account in the case of 
a tagged particle, and this was overlooked in previous calculations of heavy quark 
energy loss [19, 20]. 

Let us remark that the above energy loss process via (-u-channel) Compton 
scattering is rare but very efficient. Neglecting logarithms and writing heuristically 
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(13) 



we see that despite a large Compton mean free path oc E, the typical energy loss in 
Compton scattering is large, (iv) ~ E, and the related rate dE /dx can thus compete 
with the t-channel contribution. 



beyond leading logarithms 



In order to determine the constant beyond the leading logarithm in (11) and 
in (12), the following procedure can be used [15, 16]. Subtracting the leading 
logarithms J m2 d\t\/\t\ and J M2 d\u\/\u\ from the complete expression of dE/dx, 

the remaining integrals are dominated by either \t\ ~ m 2 D or \t\ ~ s ~ ET. An 
accurate calculation of the former integrals requires using the resummed hard 
thermal loop [17, 18] gluon propagator in the t-channel, as pictured in Figs. 3a and 
3d. The integrals dominated by \t\ ~ s ~ ET can be evaluated with bare exchanged 
propagators, but require using exact kinematics and exact expressions of the Born 
elastic scattering amplitude. Some previous attempts to calculate the heavy quark 
collisional loss beyond leading logarithm [19, 20] were misleading, for the simple 
reason that the collinear leading logarithm had been overlooked. 

Adding the leading logarithmic contributions (11), (12), and the constant c(n/) ~ 
0. 146 + 0.050 determined in [16] we obtain the heavy quark collisional loss 

/ n f \ET 2ET , J 

( 1 + f) ln ^ + 9 ln MJ + c( "4 <14) 

Note that the contribution from the collinear logarithm is reduced due to an 
associated smaller color factor [16] . 
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2.2. Untagged light parton 

Here we briefly discuss the case of a light incoming parton. When there is no 
heavy-quark tagged jet, it is impossible to know whether the detected jet arises, 
in the elastic scattering, from the scattered incoming parton of from the struck 
thermal parton. Moreover, a light incoming parton can be annihilated, due to 
the presence of partons of the same flavour in the plasma. Thus, in the untagged 
case, the loss AE in elastic scattering must be defined with respect to the leading 
parton in the final state. Defining momenta as in Fig. 3 we have (neglecting k ~ T 
compared to E) 

u<E/2 AE = P -P^ = E-E' = cu 

(15) 

u>E/2 AE = P -K' ~E-u <E/2 

Since AE < E/2, there is no full stopping with this definition. 

Thus when u > E/2, the apparent energy loss, defined with respect to the leading 
parton, is not uj as in the heavy-quark tagged case, but E — uj. Let us express the 
condition uj > E/2 using Mandelstam invariants. We have 

\t\ =2K.K' = 2k(k + Lu)(l-cos9 kk ,) 

(16) 

s ~ 2P.K ~ 2Ek{\ - cos6 kp ) 



where we used s ^> M 2 and E ^> M in the second line. When u > E/2 ^> k ~ T, 
it is easy to check from energy-momentum conservation that p and A;' are quasi- 
collinear, implying 6kk' — 6k P - Taking the ratio of the two equations in (16) we thus 
obtain 

,-, , n \t\ oj \u\ E — LO . . 

To get the last equality we used s + u + t = 2M 2 and assumed that \u\/s ~ Oil), 
which will be justified in the following. Reexpressing (15) in terms of Mandelstam 
invariants we get 

|*| < s/2 AE/E~\t\/s 

(18) 

|t| > s/2 AE/E~\u\/s 

When \t\ > s/2, the energy loss weight in (8) is simply obtained by replacing 
|*| — > \u\ in the integrand. In the untagged case, the expression (9) of the energy 
loss is thus modified to 
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We can readily see that the first term of (19) will yield the same *-channel leading 
logarithm as in the tagged case, but that the second term does not produce any 
■u-channel logarithm, due to the energy loss weight oc |w| instead of oc |*|. The term 
of \Mi\ 2 which is oc s/\u\ indeed contributes to the second term of (19) as 
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where the integral over u is dominated by |u| ~ O(s) ~ O(ET). The contribution 
(20) is subleading compared to the *-channel logarithmic contribution. 

As a result, the energy loss of an untagged parton, to logarithmic accuracy, is 
given by the first term of (14) only, 
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where Cr = Cp = 4/3 for a quark and Cr = C\ = 3 for a gluon. The expression (21) 
is Bjorken's result [4] and is specific to the case of an untagged parton. We stress the 
importance of defining correctly the observable (tagged or untagged) parton energy 
loss in order to obtain meaningful results, even to leading logarithmic accuracy. 



2.3. Running coupling 



The results (14) and (21) obtained in the fixed coupling approximation are not 
truly predictive, since the scale at which to evaluate a s is not specified. In order 



to obtain predictive results, we must repeat our calculations including the effect of 
running coupling. 

In the case of t-channel scattering (Figs. 3a and 3d), we must set the scale of 
a s to ~ C?(|£|), as dictated by renormalization in PQCD at zero temperature [21]. 
Setting this scale in the factor a 2 of the 2-^2 cross section amounts to perform a 
specific all-order resummation, corresponding physically to the final partons being 
accompanied by collinear gluon radiation. The true exclusive 2-^2 cross section 
involving charged (coloured) particles vanishes [22]. 

With running coupling, the t-channel leading logarithm in (11) becomes 

^r d 4^ r^ozm- (22) 



D 



\t\ J m 2 \t 



Using a a (|t|) = [4tt/? In (|t |/A 2 )] \ where (3 = (H - 1%)/(4tt) 2 , we find that using 
a running coupling amounts to the replacement 

a 2 In (ET/m 2 D ) — > a s (m 2 D )a s (ET) In (ET/m 2 D ) . (23) 

As noted in [21], the latter result remains finite when E — > oo. The asymptotic 
collisional loss is bounded, instead of logarithmically enhanced as suggested by the 
fixed coupling result (14). 

For the w-channel contribution (Fig. 3c), a s should be evaluated at a scale 
~ C?(|u|). Similarly to the preceding discussion, the -u-channel collinear logarithm 
appearing in (12) is modified as 

a 2 \n(ET/M 2 ) — > a s (M 2 )a s (ET) In (ET/M 2 ) . (24) 

We have already mentioned that the constant appearing next to the leading 
logarithms (in (14) or (21)) is given by integrals where the typical value of t is 
determined by either \t\ ~ m 2 D or \t\ ~ ET. Since there is no logarithmic spread 
there, the terms oc a 2 in (14) and (21) which are not logarithmically enhanced can 
be evaluated at a scale chosen arbitrarily between m 2 D and ET. 

Using (14), (21), and the modifications (23), (24) required by running coupling, 
we can write the collisional loss of a fast (E 3> M 2 /T 3> M) heavy quark as 



dE(Q) 4 . 2 , _[/ n f \. ET 2 a s (M 2 ) , ET . . 



dx 



and the collisional loss of a light energetic (E ^> T) parton as 



dE 

dx 



= C R 7ia s (m 2 D )a s (ET)T 2 

untagged 



(25) 
(26) 



We refer to [16] for a discussion of the theoretical uncertainty. 

Due to the additional logarithm in (25), when E ^> M 2 /T the collisional loss is 
(slightly) larger for a heavy than for a light (untagged) quark. Taking E — > oo and 

n f = 3, 

dE(Q)/dx 1 + A^!) <1 + A~ 115 (27) 

dE{q)/dx e^oo ^ 27 a s {m 2 D ) 27 ~ 1 ' 



3. SUMMARY 



We discussed the average collisional loss dE/dx of a heavy quark crossing a QGP, 
in the limit E ^> M 2 /T, and briefly reviewed the case of an energetic E ^> T 
light parton. For fixed a s , at leading order dE/dx oc a 2 s) with a coefficient which 
is logarithmically enhanced. A soft logarithm is present for both heavy and light 
partons. A collinear logarithm contributes additionally to the heavy quark loss, 
but not to the light parton loss, due the difference between the definition of energy 
loss for a tagged heavy quark and for an untagged light parton. Implementing the 
running of a s we obtained the corresponding results (25) and (26). 

At very high energy the heavy quark collisional loss is slightly enhanced com- 
pared to the light quark loss, see (27). However, as mentioned in section 1, for 
E ~ p T ~ (9(20GeV) <C M 2 /T, the b quark collisional (as well as radiative) loss 
should be reduced compared to the light quark case, challenging our understanding 
of heavy flavour quenching. 
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